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ABSTRACT 

This study describes suppression of flexural 
forced vibration or the self-excited vibration of a 
rotating shaft system not by passive elements but by 
active elements. Namely, the distinctive feature of 
this method Is not to dissipate the vibration energy 
but to provide the force cancelling the vibration 
displacement and the vibration velocity through the 
bearing housing In rotation. Therefore the bearings 
of this kind are appropriately named "Active Control 
Bearings . " 

A simple rotor system having one disk at the 
center of the span on flexible supports Is Investi- 
gated In this paper. The actuators of the electro- 
dynamic transducer are Inserted In the sections of 
the bearing housing. First, applying the optimal 
regulator of optimal control theory, the flexural 
vibration control of the rotating shaft and the vibra- 
tion control of support systems are performed by the 
optimal state feedback system using these actuators. 
Next, the quasl-modal control based on a modal anal- 
ysis Is applied to this rotor system. This quasl- 
modal control system Is constructed by means of 
optimal velocity feedback loops. The differences 
between optimal control and quasl-modal control are 
discussed and their merits and demerits are made clear. 
Finally, the experiments are described concerning only 
the optimal regulator method. 

INTRODUCTION 

It Is well known that a damped flexible support 
rotor system suppresses the vibrations of the rotating 
shaft system. On this point of view, the several 
studies deal with the optimum tuned condition of a 
damped flexible support system applying a flexible 
support (1,2). But It Is very difficult for this 
method to supply appropriate damping on practical 
equipment or to change a damping coefficient according 
to system variation. So there are some suggestions 
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about an active vibration control method which Is 
different from that by a passive element. G. 
Schweitzer (3) tried vibration control of a rotating 
shaft by means of an active damper using electromag- 
netic forces. In the same way, J.L. Nikolajsen et al. 
(4) carried out lateral vibration control of a marine 
transmission shaft with an electromagnetic damper. 

Host recently J. Salm et al. (5) have discussed a 
problem of vibration reduction for a rotor with mul- 
tiple degrees of freedom and have performed vibration 
control of an unbalance response using a magnetic 
bearing. H. Ulbrich et al. (6) have conducted an 
experiment which enlarges a stable region of an asym- 
metric rotor by applying a magnetic bearing. U. 
Gondhalekar et al. (7) have shown how to apply a con- 
trol force to three poles of a magnetic bearing from 
two directional signals using a microcomputer and 
have Indicated they are able to shift their critical 
speeds as a result of changing the stiffnesses. 

However the above expressed methods are not 
available when the setting positions of actuators are 
restricted. It Is complicated and dangerous to con- 
trol the shaft directly In high speed rotation. 
Therefore It Is considered to suppress the vibration 
of the rotating shaft by controlling the nonrotating 
bearing housings. R. Stanway et al. have suggested 
an active control of bearing housings and discussed 
both controllability and observability (8). Moreover 
they have studied a pole assignment problem for a 
three degree of freedom system (9). J.W. Moore et al. 
have tested the efficiency of this method by means of 
control systems using a transfer function with veloc- 
ity feedback loops (1_0) and experiments using loud 
speakers (H). 

This paper Is on the basis of an Idea similar to 
that of R. Stanway et al. Namely, a control method 
which moves the characteristic roots of a rotor system 
to the optimum positions Is applied. The efficiency 
of this control method Is demonstrated by simulations 
and experiments. As the bearing housings are actua- 
tors of the vibration control, the bearings of this 
kind are appropriately named to the "Active Control 
Bearings." In this paper two control methods for 
active control bearings are discussed In particular. 
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One 1s the control method by an optimal regulator with 
all state variable feedback ( 1 _ 2 ) and the other is by a 
quasl-modal control with velocity feedback based on a 
modal analysis. First, their design methods are 
explained about how to construct the control loop. 
Next, the two methods are compared as to their vibra- 
tion control effects, their merits and demerits are 
made clear. The rotor system has two degree of free- 
dom system with one disk at the center of the span on 
flexible supports. And the actuators were made by the 
use of a principle of an electrodynamic transducer and 
the experiments were performed by setting them. The 
results obtained agree with the simulations qualita- 
tively. Therefore it is clear that unbalance vibra- 
tions can be suppressed sufficiently by active control 
bearings. 

NOMENCLATURE 


Xr rotor absolute displacement in x direction 

Xj bearing housing absolute displacement in x 
direction 

Yf. rotor absolute displacement in y direction 

Yj bearing housing absolute displacement in y 
direction 

c rotor eccentricity 

P weighting coefficient between vibration energy 
and control energy 

Cr rotor damping ratio, fp = Cr/2-^m^ 

fs support damping ratio, = C 5 / 2 ym^ 


A system matrix 

b control vector 

Cp rotor damping coefficient 
Cj support damping coefficient 
F modal feedback gain matrix, n x n 

K stiffness ratio, K = kp/k; 

Kp stiffness matrix, n x n 

kp rotor shaft stiffness coefficient 

k; support stiffness coefficient 
M mass ratio, M = mp/mj 

Mp mass matrix, n x n 

fflp rotor mass 

mj bearing housing mass 

P forced external vector, n 

Q weighting matrix 

R solution matrix of Riccati equation, solution of 
Eq. ( 8 ) 

T modal matrix, n x n 

U control input vector, n 

Ux control input in x direction 

Uy control input in y direction 

X state vector 

Xn column vector, n 


u rotor angular velocity 

<op single support critical speed, up kp/mp 

VIBRATION CONTROL BY MEANS OF OPTIMAL REGULATOR METHOD 


Model of shaft system 

Figure 1 represents the model of a rotor shaft 
system. The shaft is considered as a massless elastic 
member and the rotor mass is concentrated in a disk 
mounted at the center of the span. The bearings are 
assumed rigid. The bearing housings are supported on 
damped flexible supports. If the control forces are 
added at both bearing housings as shown in Fig. 1, 
equations of motion for the rotor shaft system in 
complex notation reduce to the following. 


2 . 1 ut 


+ C^(X^ - Zj) - k^(Z^ - Zj) = m^e<.>‘e 

Vs ^ C3Z3 - C^(Z^ - Z^) . k^Z^ . k^(Z^ - Z^) = U 
where 


( 1 ) 


Zp = Xp + lYp, Zo 


Xs + IY 5 , u = Ux + iUy 


( 2 ) 


Considering dimensionless parameters, if only the x 
direction is shown because there is no coupling 
between x and y directions, their forms can be written 
as 


+ 2Zp“p(X^ - Xj) + co^(X^ - Xj) = J cos ut 




ZC^M^p + 2fs 


Vf"r) 


Xs - 2f^M.^X^ 


(3) 




2 ^ 
«“r^ = mJ 


where Xp/,. and Xj/^ are newly expressed as x/ 
and Xj. Transforming Eq. (3) into the state equation 
yields 


X = AX + bUj^ + wfjj (4) 

where 
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In the same manner, the equation for the y direction 
may be expressed as follows: 


Y = AY + bU + wf (5) 

^ 'S' y ~ y 

where 

Y = (Y^, Y^, Yj. Yj)"^. fy = < 0 ^ sin cot 

The superscript T denotes the transpose. The state 
Eqs. (4) and (5) mean the control systems for the case 
of one input and four outputs. Comparing two control 
inputs Ux and Uy, they are identical except for the 
phase between and Uy. Therefore the system of 

only the x direction is considered after this. 

Control system design with optimal regulator 

In this section, the control input as shown 
in Fig. 1 is determined based on the optimal regulator 
theory in the modern control fields (13). Namely for 
the state Eq. (4), the problem is to find the control 
input Ux to minimize the following cost function 







(6) 


where Q is a symmetric positive semidefinite matrix 
and it is assumed Eq. (4) is controllable. The 
optimal control input to minimize Eq. (6) is 

given with state variable feedback as follows: 


U° = - b^^ (7) 

where R satisfies the following Riccati type algebra 
equation 

RA + A^R - RW^R + Q = 0 (8) 

In Eq. (7), the following expression 

f° = b’^R (9) 

is the optimal feedback vector. Substituting Eq. (7) 
into Eq. (4), the optimal closed loop system is 
expressed except for the external force as follows: 

X = (A - bb^R)X (10) 

Equation (10) is called an optimal regulator. 


Weighting matrix and root locus 

Equation (10) is asymptotic stable in any case of 
Q and the eigen values are automatically determined 
by Eq. (8) if Q is given. Therefore it is examined 
that the eigen values of the closed loop system show 
the behaviors in what manner as the variations of the 
weighting coefficients p. Hence the parameters used 
are M = K = 1 , “r = ^00 rad/s and 

ms = 1 kg. Figure 2 represents the root loci; (a) is 
the case of 0 = diag (1,1, 1,1) and weighting coeffi- 
cients of the rotor displacement, its velocity, the 
bearing housing displacement and its velocity are all 
the same. When ^ = 0, it shows the open loop system 
or the uncontrolled system. Consequently the first 
and the second natural frequencies become the undamped 
ones and they are 0.62 and 1.62. As p increases, it 
is found that the damping ratio of the first and the 
second natural frequencies Increase and the damped 
natural frequencies decrease. In this Fig. (a), it 
shows from p = 0 to p = lO^. The root loci in the 
case which p Increases larger is considered by the 
root loci of the square of Kalman equation. From 
these results, in the case that p is infinite, three 
roots approach to equivalent zeros S = 0, S = ±1ur 
and the remaining root goes toward to the infinite 
point for the left direction along the real axis as 
like the first Butterworth pattern. 

(b) is in the case of Q = diag(1 ,1 ,0.1 ,0.1 ) . It 
draws the root loci that the weights of the rotor dis- 
placement and its velocity are ten times as much as 
the ones of the bearing housing displacement and its 
velocity. When p Increases, they separate from the 
Imaginary axis and their damping ratios Increase. But 
the second damped natural frequency is inclined to 
Increase. In the case of p = 10^, the second damped 
natural frequency increases one and half times as much 
as the one of (a) and both damping ratios are almost 
the same. When p Increases beyond this, their loci 
are similar to the ones of (a). That is, three roots 
converge to three equivalent zeros and ore root 
indicates the first Butterworth pattern, (a) in more 
sensitive than (b) concerning the sensitivity of 
behaviors of their roots. Their roots hardly move 
beyond p = 10^ in (a) and p = 10^ in (b). 


Steady state unbalance response 

In order to compute the steady state unbalance 
response extending whole rotational speeds, Eq. (3) is 
solved analytically. of Eq. (3) is expressed by 
Eq. (7) as follows; 


^x = 


fl^ - ^2\ - ^3^s - 


( 11 ) 


where f 1 .f 2 .f 3 ^nd f 4 are elements of the row 
vector of fu. Substituting Eq. (11) into Eq. (3), 
the following equation is obtained. 


M X + C„X„ + K„X„ = e cos ut (12) 
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From Eq. (12), the steady state solutions Xp are 
described as follows: — 




X|. = ^ — 2 cos(o)t - 0 .) 




(13) 


Xj = ■ " 2 ’ "g ^ cos(ut - eg) 

a + b 


where 

a = (Kv) - u 2 Hn)(K 22 - u^H 22 ) - w^Cv|C 22 - Kl 2 K 21 

+ “^Ci 2 C 2 i 

b = (Kn - a) 2 Hii)oiC 22 ( 1^22 " tj 2 H 22 )<*>Cn - I^ 12“^21 


- l< 2 io)C -|2 

P] = a(K22 “ <j 2H22) * coC22b 
?2 = <oC 223 - (1^22 ~ <>>2H22)b 
P3 = - aK2'| - <oC2ib 
P 4 = - ao)C 2 '] + bl^21 
61 = tan-T ( P2/Pi ) 

02 = tan-1 (P 4 /P 3 ) 

C,,, and K., are elements of H , C , and 
IJ IJ iJ J1 ~I1 

matrixes. Figure 3 shows the unbalance responses 

X and X^ given by Eq. (13). In this case, the 
parameters used are M = K = 1 , = 0 , 

0 )^ = 100 rad/s m^ = 1 kg and = dlag ( 1 , 1 , 1 , 1 ). 

The peak amplitudes at the first and the second crit- 
ical speeds decrease In the case when p Is large. 

If the same peak amplitudes In the case of p = lO'l 
are to be realized by passive devices, damping ratios fp 
and fs of flio™ ’’’ost b® about 0.1.' It seems that It Is 
generally difficult for rotor systems supported on 
rolling bearings to supply such large damping ratios. 

For any the second peak amplitudes are higher than 
ones of the first. This Is caused by that the damping 
ratio of the second mode Is smaller than one of the 
first mode for root loci 1n Fig. 2(a). Even If p 
Increases 100 or 1000 times as much as p = 10^, It 
can not be expected that the peak amplitudes at reso- 
nance decrease remarkably. And yet, the critical 
speeds can be actually Ignored In the neighborhood of 
p = IQb. For the case of the larger weighting coef- 
ficient p. It seems from the root loci In Fig. 2(a) 
that damping ratios of complex eigen values decrease. 
However, since the root of the second mode hardly 
moves and the unvibratlonal root on the real axis have 
an effect on the complex root, the resonance peak does 
not appear nearby the critical speed on simple sup- 
ports. This point fundamentally differs from the 
design of the optimum tuned dynamic damper and Indi- 
cates definitely that active control bearing systems 
of this kind are superior to any passive vibration 
control device. This reason Is caused by the all 
state variable feedback, specially the optimal state 
variable feedback. If an active control bearing 
system Is designed by an Incomplete state variable 
feedback, for example, only the displacement feedback 


or only the velocity feedback, the resonance peak 
arises nearby the original critical speed In 
the occasion of large loop gains. The control system 
becomes unstable In the worst case. These are dis- 
cussed In detail In next Chapter. 

VIBRATION CONTROL BY MEANS OF QUASI-MODAL CONTROL 
METHOD 

Control system design method with quasl-modal control 
First, a control system design method Is 
described In general. The following equation Is 
considered. 

Jli " J 1 Ji “ 

Equation (14) Is generally written and however a damp- 
ing term Is neglected for the simplicity. For the 
X^ vector of Eq. (14), the following transformation 

Is performed. 


^ ® = r\ C5) 

where T Is the modal matrix and It Is assumed to be 

normalized. From Eq. (15), Eq. (14) Is transformed 
as follows: 

T^M Ta + T^K Ta = T^P(t) + T^U(t) (16) 

~ „ J1™. 

Hence T^M„T Is the unit matrix and T^K T Is the 
~ n ~ n 

frequency matrix. Therefore Eq. (16) Is reduced to 
the following equation. 

a + fi^a = T^P(t) + T^U(t) (17) 

where symbolizes the frequency matrix. After 
all, Eq. (14) Is reduced to Eq. (17) separated the 
each mode and uncoupled. By the way. It Is very easy 
to determine the optimal control Inputs to each mode 
In the above expressed modal domain. That Is, It Is 
assumed to determined as follows: 

a t n^a = T^P(t) + u°(t) (18) 

where Is the new term determined as the desir- 
able optimal control force. As Eq. (17) and Eq. (18) 
should be equivalent, the following relation Is 
Introduced. 


UO = T^u (19) 

From this, the control force supplied to the actual 
system Is 


U = (T^)”%° = H"^U° 


(20) 


where ^ = T,^. Now for the undamped vibration system 
It Is c'onsidered that optimal control Inputs are 
velocity feedbacks assigned the critical damping on 
each mode. In this case, the modal control system Is 
expressed by the block diagram as shown In Fig. 4. 


where 


4 





F = 




0 -2C„w„ 

n n 


(21) 


Finally the modal control Input 1s given from Fig. 4 
by 

U = (22) 


In n degree of freedom system, If the number of 
control Inputs are less than n or the number of 
vibration velocities measured are less than n, this 
control system can not be called a modal control 
system. This paper names the such control system to 
a quasl-modal control system. For example, for two 
Inputs and three outputs, the quasl-modal control 
system 1s expressed by 
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If It Is desired that modes from the first to the s^^^ 
are controllable, matrix £ has to be regarded till 
the s^h mode. And higher ^rder modes above this can 
be disregarded. That Is, 




-2?2“2 


-2?s“s 


F = 



(24) 


Anyway the feedback gain of the quasl-modal control 
Is given by In Eq. (22) and a linear relation exists 
between U and X. 


Steady state unbalance response 

Figure 5 represents the unbalance response for 
the model In Fig. 1. The rotor vibration velocity and 
the housing velocity are measured and control forces 
are supplied the bearing housing as well as Chapter 3. 
Accordingly, this model becomes a quasl-modal control 


system with one Input and two outputs. If this system 
has two Inputs and two outputs. It Is an Ideal modal 
control system. Figure 5 shows the results for the 
several modal damping ratios f of Eq. (24). In this 
figure, displacements X^ and Xj show the mini- 
mization nearby 0.35 or 0.7. For the case that f 
Is larger than It, new peak amplitudes appear nearby 
u/<op = 0.85. The above expression leads to the 
result that the optimal damping ratio Is the vicinity 
from 0.35 to 0.7 In this quasl-modal control system. 

On the contrary, for the Ideal modal control system, 
the more the damping ratio f Increases, the more 
the amplitudes X^ and Xj decrease. Namely It has 
the most desirable characteristics of the vibration 
control. 

Comparison optimal regulator system with quasl-modal 
control system 

Figure 5 Illustrates the amplitude characteris- 
tics In the case of the optimal regulator with 
p = 104 and p = 10^. From Fig. 5, It Is obvious 
that the peak amplitude by means of the optimal regu- 
lator Is equivalent to the minimum amplitude by means 
of the quasl-modal control method. This means that 
the efficiency of the vibration control by using the 
optimal regulator method approaches to a saturation. 
Because, even If the control system Is constructed by 
the optimal control theory, it is theoretically Impos- 
sible to reduce the amplitude beyond the minimum 
amplitude of the quasl-modal control system. 

From two standpoints, a modern control theory 
and a classical control theory, the design methods of 
active control bearings have been discussed. The 
merits and demerits between the optimal regulator 
system and the quasl-modal control system are sum- 
marized as follows: 

The optimal regulator method: 
merits: (1) If a control object and a cost function 

are given, the feedback coefficients are 
automatically determined. 

(2) A closed loop system is always stable 
and optimum regions that peak amplitudes 
become minimum don't exist. 

(3) A closed loop system Is reasonable 
because the feedback control system uses 
all state variable. 

demerits: (1) It Is difficult for a multidegree of 

freedom system to give a weighting 
coefficient matrix. 

(2) It Is complicated for a raultidegree 
of freedom system to compute a feed- 
backgaln. 

(3) All state variables are required. 

The quasl-modal control method: 

merits: (1) The physical prospects are good. 

(2) It Is easy for a raultidegree of 
freedom system to determine the feedback 
coefficients . 

(3) Since this method consists of only the 
velocity feedback, the number of measure 
ments are a few. 

demerits: (1) If a damping term measurements are 

regarded. It Is complicated to make the 
modal transformation. 

(2) Optimum regions that peak amplitudes 
become minimum always exist. And for a 
large feedback gain the closed loop 
system becomes unstable. 

(3) Since It consists of only the velocity 
feedback, the control system apt to be 
unreasonable. 
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As damping forces always exist in the actual rotating 
machinery, these amounts have to be correctly esti- 
mated if possible. For the quasi-modal control system 
in particular, estimated damping amounts have sensi- 
tively an influence to the efficiency of the vibration 
control on account of the existence of the tuned con- 
dition. Occasionally the rotor system becomes 
unstable because of the miss matching of this kind. 
Judging totally from the above mentioned results, it 
is able to conclude that the optimal regulator method 
is superior to the quasi-model control method for the 
control design of an active control bearing system. 

EXPERIMENTS 

Test rig 

Figure 6 shows the overview of the test rig. The 
rotor shaft is 10 mm in the diameter, the drill rod 
with the span length 800 mm and has one disk of the 
mass 1 kg at the center of the span. The mass of the 
bearing housing is 0.5 kg including the ball bearing. 
Therefore the mass ratio between the rotor and the 
bearing housing is unit. There are four springs sus- 
pending the bearing housing and its spring constant 
is 0.49x10^^ N/m (0.5 kgf /mm) . The actuators are 
trial manufactured using a principle of a electro- 
dynamic transducer. The active control bearing set 
up by means of four actuators is indicated in Fig. 7. 
The generating force of one actuator is approximately 
18 N versus the coil current 1 A and the force is in 
proportion to the current. Figure 8 shows the sec- 
tional view of one actuator. 

Figure 9 illustrates the experimental block 
diagram. In this Chapter, it is described in only 
the case that the active control bearing system is 
designed by an optimal regulator. First the rotor 
displacement and the housing displacement are measured 
by the gap sensor. 

The first order^lag filters are used in order to 
remove some noise. X^ and Xj are obtained by dif- 
ferentiating these displacements. Next multiplying 
feedback gains fi.f 2 >f 3 >f 4 4''** adding them, the 

actuators are controlled through the servo amplifier 
supplied this state feedback signal. The control 
system in y direction is also constructed in the same 
manner . 

Discussion of results 

For M = 1, K = 0.5, Up = 100 rad/s, M; = 1 kg, 
p = lO'^ and 0 = diag (1,1, 1,1), the optimal feed- 
back gains are~glven as f^ = (-3100, 50, 3100, 

1000). Figure 10(b) shows the impulse response 
obtained by the experiments with this feedback gain. 

The efficiency of the active vibration control is well 
recognized in both the rotor and the bearing housing. 
Figure 11 shows the experimental data of the unbalance 
response with the parameter as same as in Fig. 10(b). 

In the uncontrolled case, there are two critical 
speeds at about 750 and 1600 rpm. The maximum ampli- 
tude at the first critical speed attains about 2.3 mm 
in Xf.. However, in the controlled case, the ampli- 
tude at the first critical speed reduces to about 

0.4 mm and the second critical speed can not be con- 
firmed. Figure 10(c) represents the results in the 
case of the larger weighting coefficient p, namely 
p = 3.5x10'^ and f° = (-9000, 100, 9000, 200). Under 
this case, the Impulse response of the housing is made 
better in comparison with Fig. 10(b) and the unbalance 
response is shown in Fig. 12. In spite of the 
response of the rotor in Fig. 12(a) is as same as one 
in Fig. 11(a) except for the drop of the first crit- 
ical speed, the efficiency of the vibration control 


of the housing is Improved better than Fig. 11(b). 
Figure 13 indicates the numerical solutions concerning 
the experimental conditions in Fig. 11 and Fig. 12. 

The used parameters are H = 1, K = 0.5, up = 100 rad/s 
and Ms = 1 kg. Besides the damping ratios are 
estimated as - 0.02, = 0.03 by the analysis 

of the uncontrolled experimental data. From these 
results, it can be admitted that both the experiments 
and the simulations have good agreements. 

CONCLUSIONS 

This paper proposes two control design techniques 
to the active control bearings. The superiority or 
inferiority of these design technquies has been 
described in detail by the simulations. If a rotor 
shaft system is supported by active control bearings, 
it is possible to change damping ratios of the system 
because the arbitrary state feedback system or the 
velocity feedback system can be easily constructed. 

From this, the rotor system supported by the active 
control bearings system can be designed as satisfying 
a design specification for resonance amplitudes. 
Considering the difficulty of estimation of damping 
forces in a practical rotating machinery, it is 
expected that the optimal regulator method is superior 
to the quasi-modal control method for the control 
design of the active control bearing system. The 
active control bearing system is effective not only 
for unbalance forces but for unstable forces or 
external forces transmitted from a foundation. 
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Fig. 1. - Rotor shaft model and coordinates. 
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(b) Housing displaceitient. 

Fig. 3, - Steady state unbalance response using optimal 
regulator method. 



Fig. 4. - Modal control system with velocity feedback. 
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Fig. 5. - Steady state unbalance response using quasi 
modal control method. 




Fig. 6. - Test rig overview. 


Fig. 7. - Active control bearing supported by four actuators. 
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Fig. 8. - Sectional view of one actuator. 
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Fig. 9. - Experimental block diagram . 
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(a) Uncontrolled case. 





(b) Controlled case, p - 10^. 
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Fig. 10. - Impulse response. 
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(b) Housing displacement. 

Fig. 11. - Unbalance response (experimental), p = 10^. 
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Fig. 12. - Unbalance response (experimental), p = 3.5x10^. 





Fig. 13. - Unbalance response (computed). 
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